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ABSTRACT 



We present the results of 450 and 850 //m continuum mapping of the H II region 
KR 140 using the SCUBA instrument on the JCMT. KR 140 is a small (5.7 pc diameter) 
H II region at a distance of 2.3=b0.3 kpc. Five of the six IRAS point sources near 
KR 140 were mapped in this study. Our analysis shows that two of these IRAS sources 
are embedded late B type stars lying well outside the H II region, two are a part of 
the dust shell surrounding the H II region, and one is the combined emission from 
an ensemble of smaller sources unresolved by IRAS. We have discovered a number of 
relatively cold submillimeter sources not visible in the IRAS data, ranging in size from 
0.2 to 0.7 pc and in mass from 0.5 to 130 Mq. The distribution of masses for all 
sources is well characterized by a power law A^(> M) oc M~" with a = 0.5 it 0.04, 
in agreement with the typical mass function for clumped structures of this scale in 
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molecular clouds. Several of the submillimeter sources are found at the H Il-molecular 
gas interface and have probably been formed as the result of the expansion of the H II 
region. Many of the submillimeter sources we detect are gravitationally bound and 
most of these follow a mass-size relationship expected for objects in virial equilibrium 
with non-thermal pressure support. Upon the loss of non-thermal support they could 
be sites of star formation. Along with the two B stars that we have identified as possible 
cluster members along with VES 735, we argue that five nearby highly-reddened stars 
are in a pre-main-sequence stage of evolution. 

Subject headings: H II Regions — submillimeter — stars: formation — ISM: clouds 

1. INTRODUCTION 

KR 140 (/ = 133? 425, b = 0?055; Kallas k Reich 1980) is a smah (5.7 pc diameter) H II region 
located in the Perseus arm of the Galaxy at a distance of 2.3±0.3 kpc (Kerton, Ballantyne, & 
Martin 1999). It is close to, but apparently isolated from, the W3/4/5 star formation complex. 
The region is ionized by VES 735, an 08.5 V(e) star (Kerton et al. 1999). It has been the subject 
of a recent multiwavelength study (Ballantyne, Kerton, & Martin 2000) utilizing the data set of 
the Canadian Galactic Plane Survey (CGPS; Taylor 1999). 

Surrounding the H II region are six IRAS point sources (see Table 1 and Figure 1; there is 
a faint source coincident with VES 735 as well). These are the brightest sources within several 
square degrees, suggestive of a region of enhanced star formation. In particular, infrared point 
sources located at the periphery of H II regions are of interest because they may be associated 
with star formation induced by the expansion of the H II region (e.g., Sugitani, Fukui, &; Ogura 
1991). Characterization of these sources is important in order to form a complete picture of the 
star formation history in that area. 

Determining the true nature of these sources is hampered by the low resolution of the commonly 
available IRAS image products: 4' for the Infrared Sky Survey Atlas (ISSA, Wheelock et al. 1994), 
and 1' (after HIRES reprocessing) for the IRAS Galaxy Atlas (IGA, Cao ct al. 1997) and Mid- 
Infrared Galaxy Atlas (MICA, Kerton & Martin 2000). At these resolutions a collection of objects 
can appear as a single point source. Moreover, areas of extended dust shells surrounding H II 
regions can be misidentified as point sources. While the IRAS colors can be used as a general guide 
to the nature of the source (e.g., Hughes & MacLeod 1989) they can sometimes be misleading (see 
Sections 4.1.1 and 4.1.4). 

Moving from the infrared wavelengths observed by IRAS to the submillimeter yields two advan- 
tages. First, the current generation of submillimeter mapping bolometers have far better resolution 
than the IRAS data and sufficient sensitivity to allow star formation regions at the distance of 
KR 140 to be investigated. Second, the submillimeter observations are sensitive to colder dust 
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not detected by IRAS and thus allow a more complete picture of the dust distribution. The ob- 
servations we have obtained have both clarified the nature of the IRAS point sources surrounding 
KR 140 and have led to the discovery of new discrete sources not visible in the IRAS data. 

In Section 2 the observations are described. Data analysis techniques are outlined in Section 3. 
In Section 4 the results of the analysis are discussed. Finally, in Section 5 the conclusions are 
presented. 

2. OBSERVATIONS 

The observations were obtained through the CANSERV service observing program at the 
James Clerk Maxwell Telescope (JCMT ^) using the Submillimeter Common-User Bolometer Array 
(SCUBA, Holland et al. 1999). Observations using SCUBA in its "jiggle map" mode were obtained 
in December 1997 and observations in the larger-scale "scan mapping" mode were obtained in July 
1998. 

SCUBA was used to acquire 450 jum and 850 jum data simultaneously. At these wavelengths 

the beam can be approximated as a Gaussian with a FWHM of about 8" at 450 fim and 14.5" at 
850 fim (Matthews 1999). The main beam solid angle (J^mb) can then be calculated using the 
standard formula Omb = ^-^^^^fwum (Kraus 1966). 

2.1. Jiggle Maps 

Jiggle mapping, being closely related to pointed observations, allows users to make images of 
objects on the scale of the SCUBA field of view, ~ 2'. The northwest rim of KR 140 was selected 
as the target for the initial jiggle maps because the brightest infrared source (IRAS 02160-1-6057) 
is located there (see Section 4.1.1 for the analysis of this source). Figure 2 shows the resulting 
jiggle maps (~ 2' x 2' in size). A chop throw of 2' with a position angle of 135° was used to 
obtain a reasonable sky estimate by avoiding bright infrared emission associated with the H II 
region. Beam sizes of 8.5" and 14.5" were achieved for this set of observations. The average 
normalized sky transparency during the observations was 0.81 at 850 fim and 0.34 at 450 fim. The 
expected noise equivalent flux density per pixel for the bolometer arrays is approximately 0.08 and 
0.95 Jy Hz~°-^ at 850 and 450 /um respectively (Holland et al. 1999). The measured noise levels in 
the maps are ~ 0.01 Jy/beam and ~ 0.08 Jy/beam at 850 and 450 /im respectively, consistent with 
the expectations from sky transparency and bolometer characteristics. There are no strong point 
sources visible in either map. Instead, extended emission is quite evident in the 850 //m map (S/N 
~ 4-8). The same region of emission is barely visible in the 450 /xm map. However, convolving the 



^The JCMT is operated by the Joint Astronomy Centre on behalf of PPARC of the UK, the Netherlands OSR, 
and NRC Canada. 



-4- 



image to the 850 /xm resolution reduces the noise level to 0.05 Jy/(850 /xm beam) and results in 
the extended emission being clearly visible (S/N ~ 3-7). 

2.2. Scan Maps 

Larger scale maps of KR 140 at 450/850 ^m were obtained in July 1998 using the recently (at 
the time) developed "scan mapping" observing mode which allows larger areas to be mapped out 
by scanning the array across the sky. These maps included the jiggle- map region and enabled the 
study of the other IRAS point sources near KR 140. 

The raw bolometer data were reduced using the standard SCUBA software (Holland et al. 
1999) for flatfield and extinction corrections. The data were then further reduced and transformed 
into the final maps using a matrix inversion technique (Johnstone et al. 2000a) which effectively 
minimizes the difference between the individual chop measurements and the resulting map. One 
advantage of this technique compared to standard reductions is that the bolometer input may be 
weighted to account for variations in the quality of the measurement. Figure 3 shows maps created 
at this stage. Weiner filtering these maps (see, e.g., Press et al. 1992) removes the high (spatial) 
frequency noise associated with the individual chop measurements without affecting structure on 
scales equal to or larger than the JCMT beam. Further processing of the map proves useful. Due 
to the finite chop size, SCUBA maps have their low spatial frequency modes suppressed. Restoring 
these modes requires a large re-amplification which also amplifies any existing low frequency noise. 
Furthermore, correlated errors introduced by variations in the sky opacity and the drift in the 
bolometer response produce excess residual low frequency modes. Thus the large scale structure in 
the map is less reliable, and in order to better visualize the small scale features in the data, the low 
frequency modes were entirely suppressed. This was achieved by convolving the filtered map with 
a cr = 130" Gaussian beam (twice the largest chop throw) and then subtracting the convolved map 
from the original to remove any large-scale structure. Both images were then brought to a common 
20" resolution to facilitate comparison. In Figure 4 the resultant maps are displayed. Noise levels 
are ~ 0.01 Jy/beam and 0.05 Jy/beam at 850 and 450 /um respectively. 

Figure 5 compares the derived 850 fim scan map for the region with other data obtained from 
the CGPS database (see Ballantyne et al. (2000) for details). 

3. ANALYSIS 

Our quantitative analysis focuses on the discrete sources detected in our maps. The next three 
subsections describe how the sources were detected, and the subsequent measurements of flux, size, 
temperature, mass and density for these objects. Interpretation of these results is presented in 
Section 4. 
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3.1. Source Detection, Flux and Size 

Individual sources were identified through a combination of visual inspection of the maps, and 
by using a 2D version of the automated clump finding algorithm CLUMPFIND (Williams, de Geus, 
k, Blitz 1994) which separates objects along minimum flux surfaces within the map. In total 22 
sources were detected; 11 at both 850 and 450 /im, nine sources at only 850 /im, and two at only 
450 ^m. These last two sources may be visible at 850 //m but it is hard to be certain because of 
the increased noise at that location due to edge effects from processing. The flux-weighted centroid 
of each source was determined and used as a measure of the source position. We have assigned 
the acronym KMJB to these new submillimeter sources. Except for one of our sources (KMJB 1 
= IRAS 02174+6052) detected near the edge of the map, there was no significant difference in the 
centroid positions between the two observed wavelengths (when applicable). For KMJB 1 we have 
reported the 850 ^m centroid position because of the higher quality of the data. Table 2 lists the 
positions of the sources indicated by the white crosses in Figure 4. 

In order to estimate the size of the sources, an azimuthally averaged Gaussian profile was 
fitted to each source using the image analysis program kview (Gooch 1995). To take into account 
the beam size, a simple convolution of a Gaussian source with a Gaussian beam was assumed, 
^observed = ^beam + ^object' ^^^re 6'6eam = 20". Physical sizes {D^s = Oohiectd) were calculated 
assuming a distance d = 2.3+0.3 kpc to KR 140 and tabulated in Table 2. The effective diameters 
range from 0.7 pc down to the resolution limit, 0.1 - 0.2 pc. This size scale of a few x 0.1 pc is 
typical for cores observed within molecular clouds (Evans 1999). 

Flux measurements of the sources were obtained using the Imview software package^. The 
uncertainty in each measurement is dominated by the uncertainty in choosing a correct background 
level. To quantify this uncertainty, we interactively selected three different bounding polygons 
for each source, and used two different functional forms for the background (twisted plane and 
twisted quadratic) which were fit to the boundary points of the polygon. The six resulting flux 
measurements were then averaged together and the standard deviation of the measurements was 
calculated and used as a conservative uncertainty estimate. The results are summarized in Table 2. 

3.2. Temperature 

The spectral energy distribution (SED) of the dust emission was parameterized using what we 
call a "blucbody" spectrum after Dent, Matthews, & Ward-Thompson (1998): a single temperature 
Planck function with a v'^ dependence (/3 = 2) for the dust cmissivity While this is clearly a 
gross oversimplification of the dust properties Dent et al. (1998) and Gordon (1995) point out 



^Imview is a multi-purpose image display and analysis package developed at DRAO. 

^ "Modified blackbody" is the term often used in the literature, but bluebody is both shorter and more descriptive. 
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that the quantities derived from these simple fits are often consistent with results derived from a 
more complex, and often non-unique, analysis of the data. 

For the isolated sources detected by IRAS, this SED was fit to the submillimeter and 100 and 60 
/xm IRAS data. However, for two sources, IRAS 02168+6052 and IRAS 02157+6053, this technique 
was not possible because of the significant beam size difference between IRAS and SCUBA: multiple 
submillimeter sources were detected within the larger IRAS beam. In these cases, and for sources 
without IRAS counterparts, the submillimeter photometry was used where possible to derive a 
450/850 /xm bluebody color temperature. This estimate can be rough for two reasons. For warmer 
dust, one is out on the Ray leigh- Jeans tail of the Planck function which is relatively insensitive to 
temperature. The calculated temperature for colder dust is still very uncertain due to the relatively 
poor quality of the 450 /xm data. 

For sources detected at only a single wavelength representative masses were calculated using 
a temperature of 20 K, as has been done in other single wavelength analyses of millimeter and 
submillimeter data (e.g., Johnstone et al. 2000b, Motte, Andre, & Neri 1998). This temperature 
is a typical value observed for compact dust features found by submillimeter studies of the Orion 
molecular clouds (Chini et al. 1997; Lis et al. 1998). 



3.3. Mass and Density 

Table 2 summarizes the mass and density estimates. The mass calculation used the standard 
formulation of Hildebrand (1983) for an optically thin molecular cloud with a uniform temperature: 

M = d^F^C^/B^T) (1) 

where M is the total (dust and gas) cloud mass, d is the distance to the objects, Fi, is the observed 
flux density, By{T) is the Planck function evaluated at dust temperature T, and is a factor 
combining the dust opacity {ny) and the dust to gas mass ratio {Cy = Mg/(MdKi,)). 

Table 3 summarizes various values of Cgso derived from the literature. Most of the articles cited 
report the value of at 1 mm, and the values shown have been calculated assuming that oc 
{(3 = 2 for the opacity). The canonical Hildebrand (1983) value is derived from the observations 
of the reflection nebula NGC 7023 (Hildebrand 1983). While this is empirically derived, ideally 
it should be applied to observed objects which have similar properties to the calibration object 
(Hildebrand 1983); this is probably not the case here where we find lower temperatures and higher 
densities. The Draine & Lee (1984) value is appropriate for dust in the general ISM; however, it is 
known that properties of dust in dense clouds differ from those in the diffuse ISM (Kim, Martin, & 
Hendry 1994). The difference in optical properties arises from a combination of coagulation, leading 
to an increase in grain porosity, and the accretion of volatile grain mantles. The remaining values 
in the table are from various models that have attempted to deal directly with these processes. 
Various authors suggest different values ranging at 850 /xm from 21 - 86 g cm~^, depending upon 
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the gas density and temperature. Since we do not know the density in advance we have adopted 
a value of = 50 g cm~^ at 850 /im and, assuming a (3 = 2 dust emissivity, = 14 g cm^^ 
at 450 ^m for this investigation. This is in general agreement with the results of a recent detailed 
observational study of circumstellar material by van der Tak et al. (1999). 

The masses derived in our study range from 0.7 to 130 Mq. They can be linearly scaled to 
another value of if desired. The uncertainty in the distance to KR 140 leads to a systematic 
26% uncertainty in the derived masses. The average tabulated uncertainty in the flux density mea- 
surements is 17% at 850 /xm and 28% at 450 /xm. Uncertainties relating to the derived temperature 
are harder to quantify because of the simple bluebody model that we have used and because the 
temperature enters into the mass derivation non-linearly through the Planck function. Figure 6 
shows the change in the mass derived from 850 fim flux measurements produced by varying the 
temperature. At lower temperatures the derived mass is very sensitive to the temperature. For 
example, reducing the temperature by only 2 K will cause the derived mass to almost double. The 
same absolute change at 15 K causes an increase of only 28% and at 25 K the change is only 12%. 
Based upon these considerations a total uncertainty (not including the uncertainties in Ci,) for the 
derived masses of ±50% is an optimistic minimum value (see, e.g., Gordon 1995). 

For each source the cube of D^s was taken as a rough estimate of the volume and used to 
calculate a density. Considering both the lack of information about the true three-dimensional 
structure of the sources and the uncertainty in the masses, these densities should be considered as 
only order of magnitude estimates. 



4. DISCUSSION 

In this section we first consider what the new submillimeter observations can tell us about the 
nature of the IRAS point sources, drawing on the multifrequency data available. Discussion about 
the newly detected sources and the region as a whole follow in Sections 4.2, 4.3 and 4.4. 



4.1. The Nature of the IRAS Sources 

4.1.1. IRAS 02160+6057 = KMJB 15 

The source IRAS 02160+6057 identified with the north-west dust arc has been the subject of 
two molecular line investigations. Wouterloot & Brand (1989) identified it as a potential star- 
forming area via its IRAS colors (see their paper for the exact selection criteria) and examined 
it (along with about 1300 other IRAS point sources) for CO emission. They found a CO feature 
in that direction at a velocity of —49.7 km (LSR), which corresponds to CO in the associ- 
ated background molecular cloud (Ballantyne et al. 2000), and assigned it the catalog number 
[WB89]417. This point source was not detected in a subsequent H2O line survey by Wouterloot, 
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Brand, & Fiegle (1993). 

A comparison of the scan maps with the IRAS images of KR 140 for this region provides a 
striking example of how much the morphology of the dust emission associated with an H II region 
can change with wavelength. In the IRAS images this region at all bands (Figure 5; Ballantyne 
et al. 2000, Figure 2) is the dominant feature. In contrast, while the feature is visible in the 
submillimeter maps it is certainly not the dominant one. This contrast indicates that the dust 
emission one sees from this area is from relatively hot dust with no additional cool components 
along the line of sight. This is consistent with the lack of significant clumped CO emission in 
this area which could have indicated a region containing colder dust grains (see Figure 5). The 
submillimeter and infrared emission are also aligned spatially which is what we would expect if the 
submillimeter emission is just coming from the hot dust we see in the IRAS bands (see Figure 2). 
Morphologically, IRAS 02160+6057 appears to be part of the partial dust shell surrounding KR 140. 
This line of sight has one of the highest column densities of warm dust in the KR 140 H II region 
and so it is possible that a protostar could be forming there as a result of the expansion of the H II 
region. However, examination of the HIRES images shows no point-like features in the dust arc 
and none are found in the new submillimeter observations. 

The average mean flux density at 850 /xm is 0.7±0.1 Jy in the jiggle map and 1.0±0.3 Jy in the 
scan map. For our calculations we have used 0.8+0.2 Jy as the flux density at 850 fxm. At 450 fim 
we used the flux density from the jiggle maps, 4 ± 1 Jy, since the feature is not seen in the 450 /xm 
scan maps. In order to construct an SED for the source we measured infrared flux densities in the 
same region at 60 and 100 /xm using HIRES images of KR 140. We found Fgo = 42 it 5 Jy and 
Fioo = 101 ± 8 Jy. Combining this infrared data with the submillimeter data we obtain a best-fit 
bluebody curve to the SED with T= 27.5 K (see Figure 7). Comparison with the azimuthally 
averaged dust temperatures for KR 140 (see Figure 5 in Ballantyne et al. 2000), shows that our 
derived temperature is consistent with dust being heated by the central star VES 735. 

We conclude that IRAS 02160+6057 is not a point source but simply part of the dust shell 
associated with KR 140. This naturally explains the non-detection of H2O in the Wouterloot et al. 
(1993) study. Using equation 1 we find a mass of 9 Mq in this segment of the shell; the density 
there is 10^-^ cm~^, the lowest for all of the sources and would be even lower if allowance were made 
for limb brightening (line of sight distance > -Deflf)- The structure of this dust shell (morphology, 
density, temperature) was examined in detail in Section 6 of Ballantyne et al. (2000). 

4.1.2. IRAS 02168+6052 = KMJB 6, 7 and 8 

At submillimeter wavelengths the region around the single IRAS point source IRAS 02168+6052 
contains three sources (see Figure 8). The brightest submillimeter source in the area (KMJB 8) 
is well off the 60 ^va. centroid and is associated with the CO emission. It is also detectable at 
450 //m. The resulting color temperature is 8+2 K. This leads to a mass of 80 M© and a density 
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of 10^-^ cin-=^. 

KMJB 6 breaks up into three discrete objects at 450 /xm. We calculated a 450/850 color 
temperature using the combined flux of the sources at 450 and the measured flux density at 850 
fim. The derived dust temperature is 7ib2 K. The corresponding mass is about 55 Mq, and the 
overall density in this region is around IC^'^ cm"'^. KMJB 6 and 8 are most likely cold molecular 
cores; however, the very cold temperatures derived for these objects should be confirmed with 
future spectral line observations. 

KMJB 7 is located near the position of IRAS 02168+6052 and is most likely a dense knot of 
dust and gas within the dust shell surrounding KR 140. Like in KMJB 15, the shell is the source 
of the prominent ridge of emission seen in all of the IRAS emission bands. Faint submillimeter 
emission from the dust shell is visible as a diffuse ridge running between KMJB 6 and 8. The 
temperature of dust within this shell was measured to be around 30 K (Ballantyne et al. 2000), 
which is consistent with the ratio of 60 and 100 fim to 850 /xm fluxes. Using this dust temperature 
we derive a mass of 1 Mq for the knot that is KMJB 7. 

4.1.3. IRAS 02157+6053 = KMJB 17, 18 and 19 

This region (see Figure 9) is associated with molecular material located behind KR 140 (Bal- 
lantyne et al. 2000). The three submillimeter sources are located within the larger CO envelope 
running roughly SW-NE, following the orientation of the 60 /xm contours. The IRAS scan pattern 
for this region places the major axis of the detectors along a SW-NE direction and so the PSC 
fluxes are the contribution of the ensemble of sources in the area. The crowded nature of the region 
makes it difficult to use the IRAS fiux densities to determine temperatures. 

Nevertheless, KMJB 19 is well off the 60 /xm centroid and hence relatively cold. It is detectable 
at 450 /7,m, leading to a 450/850 fim color temperature of 8ib2 K. The resulting mass (50 Mq) and 
density (10^'^ cm""^) are consistent with the object being a cold molecular core, like KMJB 8. 

Mass estimates were obtained for KMJB 17 and 18 using the 850 fim flux and assuming a 
temperature of 20 K (which is roughly consistent with explaining the 60 /xm emission). The mass 
and lower limits on the density (the sources are point sources) are both consistent with the objects 
being molecular cores. KMJB 18 is probably warmer than KMJB 17 given the facts that they are 
of comparable brightness at 850 /tm and 18 is located closer to the centroid of the 60 /tm emission. 
Taking this into account would accentuate the mass and density differences. 

4.1.4. IRAS 02171+6058 = KMJB 3 

This source is located to the northeast and outside of the KR 140 radio continuum emission. It 
seems likely that this relatively bright source is associated with molecular material at the distance 
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of KR 140. On the basis of the SED, this would be a class I object (Andre, Ward-Thompson, & 
Barsony 1993). We proceed with the analysis assuming a distance of 2.3 kpc and show that this 
consistently explains a number of observational characteristics of the source. 

The infrared colors of the infrared point source IRAS 02171+6058 are characteristic of an ultra- 
compact H II (UCHII) region based upon the criteria of Wood & ChurchwcU (1989): log{Fiy{()0)/F,y{12)) > 
1.3 and log(F^,(25)/F^(12)) > 0.57. They also match the criteria of Hughes & MacLeod (1989) for 
an H II region: log(F^(60)/F^(25)) > 0.23 and log(F^(25)/F^(12)) > 0.40. It should be noted that 
the PSC 60 fiia flux is only of moderate quality (as defined in the IRAS Explanatory Supplement 
1988), but it is in agreement with a flux density measurement we made from the HIRES images of 
the region. The source does not meet the additional flux criteria set by Hughes & MacLeod (1989), 
Fioo ^ 80 Jy which was designed to help filter out reflection nebulae and extragalactic objects. 
The source was included in a CS (2-1) survey by Bronfman, Nyman, &; May (1996) and a methanol 
maser survey by Lyder & Gait (1997), but in both cases no detections were made. 

Since this IRAS point source is detected at 12, 25, 60, 100 and 850 /im it is surprising that it 
is not seen at 450 /zm. We fit a single-temperature bluebody to the 850, 100 and 60 /xm data points 
to obtain a temperature estimate of 27 K (see Figure 7). This best fit predicts a flux density of 2.8 
Jy at 450 /xm which normally should be detectable, especially in the smoothed scan map. However, 

there is a large gradient in the 450 fim. scan map around this position due to enhanced noise at the 
edge of the map which may be making it hard to detect the source. The fitted SED yields a total 
mass of 5 Mq. Using 0.23 pc as a typical length scale, the density is found to be lO^'"*^ cm~^. 

Note that this dust is not heated by the exciting star of KR 140; there must be a source 
internal to the cloud. The total flux derived from the best-fit SED gives a luminosity of L = 10^'^ 
L0 assuming a distance of 2.3 kpc. This can be compared to the bolometric luminosity for main 
sequence B5 V - B8 V stars which ranges from L ~ 10^ ° to 10^-^ L© (Lang 1992). 

This is consistent with the lack of a corresponding point source in the 1420 MHz continuum 
image. The sensitivity of the 1420 MHz CGPS images is 0.2 mJy/1' beam (Taylor 1999). Ionization 
equilibrium can be used to relate this brightness level to the rate at which a star is emitting ionizing 
photons (Q). At the distance of KR 140 this brightness level corresponds to log(Q) = 44.0 for a 
point source, and so any embedded O or early B main sequence stars (which have log(Q) > 47; 
Kerton 2000) would certainly be detected in the CGPS radio data. 

Chini, Kriigcl, &: Wargau (1987) derived a relationship between the total luminosity (as seen 
reradiated by dust) of a compact H II region and the gas mass: 

L = (56 ± 21)M°-93±°-°6 (2) 

where L and Mg are the luminosity and gas mass in solar units. The lowest luminosity and mass 
considered in their study was L = 10^'^ Lq. Using an extrapolation of this relationship, the 
luminosity derived for this object above corresponds to a mass of 6^2 Mq. This is consistent with 
the independently determined measurement of the mass based upon the submillimeter observations. 
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and indicates that the Chini et al. (1987) relationship may hold to lower L and M values than 
covered in their study of regions with ionizing stars. 

4.1.5. IRAS 02174 +6052 = KMJB 1 

The IRAS PSC fluxes for this object are only of marginal use since they contain two upper 
limits at 60 and 100 ^m. To get more information we examined HIRES images of the KR 140 
region. At 100 (im. the source is not visible against the very strong emission coming from the 
nearby dust shell of the H II region. We were able to measure flux densities at the other three 
IRAS bands. At 60 /xm we obtained a flux density of 12ib2 Jy. This is consistent with the upper 
limit reported in the PSC. Flux densities obtained at 12 and 25 /im agreed with the PSC values. 
Like KMJB 3, this has a class I SED. The object docs have F^(60)/Fj.(25) and F^(25)/F^(12) 
colors in the range expected for an UCHII region as set by Hughes & MacLeod (1989), although 
the Fy{2b) / Fy{l2) color is at the lower limit. Wood & Churchwell (1989) used the Fy{m)/F^{l2) 
color and a higher limit to the Fi^(25)/F,^(12) color to define UCHII regions, and this object does 
not match their criteria for an UCHII region. Therefore, from the IRAS colors alone it is unlikely 
that the object is an UCHII region. 

We fit a single-temperature bluebody to the 850, 450 and 60 jiva data points to obtain a 
temperature estimate of 22 K (see Figure 7). The best-fit curve predicts a flux density of 68 Jy at 
100 /im which is consistent with the upper limits reported in the PSC. This value is also below the 
100 /xm flux density limit set by Hughes & MacLeod (1989). Values from this best-fitting curve 
were then used to derive a total mass of 19 Mq. 

Using the SED, and assuming a distance of 2.3 kpc, we calculated a total luminosity of 10^'® 
Lq. As with KMJB 3 there is no corresponding 1420 MHz continuum source. This is what one 
would expect for a star with luminosity corresponding to B5 V. Using this luminosity and the 
mass-luminosity relation of Chini et al. (1987) (equation 2) we obtain a mass estimate of 14^5^ 
Mq in agreement with the value derived from the submillimeter observations. 

4.2. Cumulative Clump Mass Function 

The cumulative mass distribution for the sources detected in this study is shown in Figure 10. 
A least-squares fit to the data of a A^(> M) = NoM~" power law resulted in a = 0.49 ± 0.04. 
Kramer et al. (1998) compiled results from a number of studies and showed that a power law 
distribution with a = 0.65 it 0.10 describes the clump mass distribution in molecular clouds for 
clumps with masses ranging from 10~^ to several 10^ Mq. The two studies contained within Kramer 
et al. (1998) that are closest to our study with respect to the mass range covered agree within the 
stated uncertainties. Clumps in the mass range 0.6 to 160 Mq surrounding the H II region S 140 
were found to follow an a = 0.65 ± 0.18 power law, while clumps in the mass range 0.8 to 50 Mq 
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surrounding the H II region NGC 1499 followed a = 0.59 ± 0.18. In Figure 11 we have plotted the 
cumulative mass distribution for the larger scale cores (0.05 pc to 0.3 pc in size) in p Oph using 
the data from Motte et al. (1998). These objects also follow a power law relationship similar to 
what wc find in KR 140, a = 0.6 ± 0.08. 

In contrast, millimeter and submillimeter studies of the nearby p Oph region have shown that 
the mass distribution of the smaller scale fragments (2400 to 5000 AU in size) contained within the 
cores appears to mimic the stellar IMF (Motte et al. 1998; Johnstone et al. 2000b). Testi k. Sargent 
(1998) used the OVRO interferometer to study similarly sized objects in the more distant Serpens 
molecular core. They also found a steep power law distribution consistent with the stellar IMF, 
a = 1.1. The key similarities in these three studies are that they probed sufficiently small linear 
scales (0.01 - 0.02 pc) and that a sufficient range in masses was observed to detect a steepening of 
the mass function above 0.5 Mq. These both appear to be important requirements. Of the studies 
listed by Kramer et al. (1998) only the study of the nearby {d = 65 pc) quiescent cloud L1457 has 
sufficiently high angular resolution to be comparable to the p Oph and Serpens studies. However, 
the highest mass clump detected is only 0.2 Mq and so all of the detected objects lie on a shallow 
part of the mass distribution function (a = 0.77 it 0.3). One can interpret this as either being 
consistent with the low mass portion of the stellar IMF or being consistent with a global power 
law for molecular cloud structure. Our study of KR 140 has sufficient overlap in mass coverage to 
probe the connection between cloud structure and the stellar IMF, but the larger scale objects that 
we detect at our relatively low resolution follow the well-established a ~ 0.5 power law seen for 
structures within giant molecular clouds and for the clouds themselves (Williams, Blitz, Sz McKee 
2000). 



4.3. Nature of the Sources 

How might the submillimeter sources relate to future star formation? The Bonnor-Ebert (BE) 
critical mass (Ebert 1955; Bonnor 1956) is given by 

2.4i?a2 

Mbe = (3) 
Using R = Des/^ and the standard formula for the isothermal sound speed (og) we can write 

M.. = 0.9.(^)(%^). (4) 

To derive Mbe from our data we have assumed that Tgas = T^ust since our CO data are of too 
low resolution to give information about the gas dynamics or temperatures within specific clumps 
directly. 

Figure 12 plots the derived submillimeter mass M versus the corresponding Mbe for each 
clump. Clumps with M < Mbe are not gravitationally bound or unstable. Prominent among 
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these are the sources that we have identified as parts of the dust shell, 7 and 15; their existence 
depends on pressure of the H II region and its dynamical expansion. All of the other "stable" 
sources (2, 10, 13, 16, 20, 21, 22) would also be pressure supported or transient. There arc clumps 
which are detected at only one submillimeter wavelength and have no clear IRAS counterparts. For 
these we have assumed T = 20 K. Adopting a lower T would have led to a lower M^e and higher 
M (Figure 6), moving the clumps toward the region of instability. However, unless we adopted a 
substantially lower temperature (~ 7 K) the clumps would remain in the stable region. 

Clumps with M > Mbe are massive enough to collapse under the effect of self-gravity. Ex- 
cluding the clumps that have already formed stars (KMJB 1 and 3), we find eleven, including the 
aligned sources 11, 12, and 14, and 4, 5, 6, 8, 9, and 19. This result derives from the low temper- 
ature (T < 13 K) that we have measured for these sources. KMJB 17 and 18 are also barely over 
the instability line. For these we have assumed T = 20 K. Adopting a lower T would move the 
clumps further into the region of instability. 

By combining the available observational data on molecular clouds, Larson (1981) showed 
that over a wide range in sizes gravitationally-bound molecular clumps follow the relationships 
a DC D^-^ and a oc M°-^ were a is the velocity dispersion, D is the largest linear dimension (in 
pc), and M is the mass of the object; thus M (x D^. This is equivalent to saying that the objects 
have a constant column density. Figure 13 shows M/D'^ (or equivalently Ay) vs. D. The bimodal 
distribution of the sources is apparent. All of the unbound sources, with M < Mbe, have a very 
low column density below Ay ~ 2. The average column density, nD = 2.1 x 10^^ cm~^, for all of 
the sources with M > Mbe (except the outlying sources 5 and 12) is twice what would be deduced 
from Larson (1981). This may be due to our using a different tracer, and thus a different length 
and mass scale. Myers (1983) showed that this result can be interpreted as virial equilibrium with 
non-thermal pressure support from turbulent gas motions. 

As discussed by Evans (1999), star formation appears to occur only in cores with Nu > 8 x 10^^ 
cm~^ {Ay ~ 4). Theoretically this is expected because of the sharp decrease in the ionization 
fraction within such cores which in turn allows ambipolar diffusion to proceed at a faster rate 
(McKee 1989, 1999). In clouds where both magnetic and turbulent support are important one still 
needs a low ionization fraction to allow collapse to proceed after the loss of turbulent support (Shu 
et al. 1999). The sharp decrease in the ionization fraction at Ay ~ 4 may also be related directly to 
the loss of turbulent magnetohydrodynamic support because of the reduction in ion-neutral coupling 
(Ruffle ct al. 1998). Onishi et al. (1998) studied 40 molecular cloud cores in Taurus in C^^O and 
found that all of the cores with A^h > 8 x 10^"*^ cm~^ were associated with H^'^CO^ emission (which 
they interpret as being protostellar condensations) and/or cold IRAS sources. In our sample of 
sources around KR 140 only KMJB 1 and 3 have unambiguous evidence for star formation having 
occurred and both are high column density sources. It would be interesting to investigate the 
other observed sources for evidence of star formation via higher resolution observations of dense 
gas tracers. KMJB 12, lying well above the line, is an interesting object for further study because 
of its large mass and location on the ridge boundary discussed below. Among the clumps with 
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M > Mbe, KMJB 5 has an unusually low column density. 

4.4. Clump and Star Formation in the KR 140 Region 

Are the clumps just structure in the cloud that pre-existed the formation of VES 735 and 

KR 140 or is there ongoing evolution? The three high density sources 11, 12, and 14 are of 
particular interest. Their cold temperature, small size, and high density suggest that they are 
cold molecular cores. They are aligned along the edge of a CO cloud with Vlsr = —46.0 km s~^ 
(Figure 5) and also parallel to, but outside of, a ridge which appears in a high resolution Ha map 
(Kerton et al. 1999). These sources are probably not foreground to the H II region since even 
with their large column densities they do not produce distinctive silhouettes on the Ha map. This 
and the striking alignment suggest that a causal relationship exists between the formation of these 
clumps at this location and the expansion of the KR 140 H II region. We have estimated that the 
age of the KR 140 H II region is about 2 X 10^ years (Ballantyne et al. 2000), which is certainly 
sufficient time for a layer of gas to be swept up and become gravitationally unstable even for the 
low density molecular cloud associated with KR 140. Other cold submillimeter sources (e.g., 6, 8, 
and 19) lie at the the edge of the H II region (see Figure 5) and their location is also suggestive of 
a link between the expansion of the H II region and their formation. 

This mix of clumps, some that appear to have formed in particular locations because of the 
action of other stars, some that appear to have formed spontaneously in various regions throughout a 
larger cloud, and some which are forming stars (this would include the clump that formed VES 735), 
is not unlike what is seen in other star- forming regions (e.g., for p Oph, see the discussion in Section 
6.3 of Motte et al. 1998). 

It would be expected that a cluster would be forming near VES 735 and so it is interesting 
to examine the stellar (or proto-stellar) content. Our analysis has shown that there are no UCHH 
regions in the vicinity of KR 140 and so VES 735 appears to be the only O star to have formed in 
the region. From Ballantyne et al. (2000) the total mass of the molecular cloud associated with 

KR 140 is estimated to be about 5x 10^ Mq. O star formation in clouds with a cloud mass of below 
10^ Mq is considered rare (Elmegreen 1985; Williams & McKee 1997) and so it would actually 
have been very unusual to find yet another O star in the vicinity of VES 735. 

KMJB 1 and 3 contain embedded objects with the luminosity of main sequence B stars. Their 
position in the molecular cloud suggests no causal relationship with the H II region and instead 
they could have formed in the original molecular cloud along with VES 735. The evolution time 
for a 4-6 Mq star (~ B5 V to B8 V) to the main sequence ranges from 0.4 - 8.4 x 10^ years (Palla 
& Stahler 1993). Ballantyne et al. (2000) estimate that VES 735 and KR 140 are 1-2 xlO^ years 
old, ample time for these stars to have formed and evolved to their present state; in fact their 
evolution is probably relatively delayed rather than coeval. 

Close to KR 140 there are five bright infrared stars (BIRS 128 to 132; Elmegreen 1980; see 
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Figure 14 in Ballantyne et al. 2000) with very red ruR — mj colors (4.3 to 6.7; m and M in this 
paragraph refer to apparent and absohite magnitudes respectively). BIRS 128 is just outside of 
KMJB 7 (see Figure 8), BIRS 131 and 132 are near KMJB 16, and BIRS 129 and 130 are seen 
projected against the H II region, north of VES 735. However, while Elmegreen (1980) considers the 
stars as possibly "embedded" there is no significant submillimeter (or infrared) emission coincident 
with their positions, which puts a limit Ay ^ 1 or Er^j ^ 0.27 on associated clumped dust. The 
general foreground extinction seen to affect VES 735 and KR 140 amounts to Ay ~ 6 (Kerton et 
al. 1999) or Er-i ~ 1.6. They are possibly more highly extinguished by the molecular material 
that lies behind KR 140 {V = —46 to —54 km s~^). Using the CO and H I data cubes and gas 
to dust conversion factors from Digel et al. (1996) or Strong &: Mattox (1996) Er^i ~ 3 seems 
possible. Nevertheless, the stars probably have intrinsically redder colors than early main sequence 
stars {Mr — Mj ~ 0). They are not ionizing stars. They could be pre-main sequence stars or 
evolved giants, though the latter seems less likely since these stars would have had lower mass 
than VES 735 yet have evolved sooner. Additional photometry at J, H, and K, when available 
from 2MASS, should help to clarify both this and the overall stellar content of the recently formed 
cluster. 



5. CONCLUSIONS 

Based on the analysis of the SCUBA 450 and 850 /xm maps of the KR 140 region we conclude 

that: 

1) The putative UCHII region IRAS 02160+6057 is actually part of the partial dust shell 
surrounding KR 140, to the NW. IRAS 02168+6052 is also part of the partial dust shell, to the 
east. 

2) The infrared point source IRAS 02157+6053 appears in the submillimeter as a number of 
discrete sources or clumps, some colder than the IRAS source. 

3) The cumulative mass distribution of clumps follows a N{> M) = NoM~" curve with 
a = 0.49 ± 0.04. This is consistent with studies of clumps of this scale in other molecular clouds. 
At the distance of KR 140 we lack sufficient resolution to study the mass distribution of even 
smaller cores, STich as has been done for the more nearby star forming regions of p Oph (Motte 
et al. 1998; Johnstone et al. 2000b), and Serpens (Testi & Sargent 1998). The use of current 
(OVRO) and future (SMA, ALMA) mm arrays is needed to build up a significant sample of star 
forming regions and to probe different star forming environments at sufficient resolution to probe 
the possible relationship between observed structures in molecular clouds and the observed stellar 
IMF. 

4) The evolution of some of the clumps appears to have been affected by the KR 140 H II region. 
Three molecular cloud cores are aligned along the edge of a larger CO structure at the interface 
between the molecular and ionized gas. The alignment strongly suggests a causal relationship to 
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the expansion of the H II region. 

5) Many of the the submilhrneter sources we detect are gravitationally bound. Most of the 
sources with M > Mbe follow a mass-size relationship expected for objects in virial equilibrium 
with non-thermal pressure support. Upon the loss of non-thermal support they could be sites of 
star formation. 

6) A cluster of stars is probably forming along with VES 735. IRAS 02171+6058, previously 
thought to be an UCHII region on the basis of its PSC colors, has the luminosity of a lower mass 
embedded B star (later than B5 V). IRAS 02174+6052 is also an embedded star with similar 
luminosity (~ B5 V). Both objects appear to have formed spontaneously in the molecular material 
outside KR 140. Their evolution is somewhat delayed relative to the evolution of the O star 
VES 735 and its H II region. We show that five BIRS stars are possibly pre-main sequence stars. 
More infrared photometry, such as will be available from 2MASS, is required to delineate the cluster 
properties. 
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Fig. 1. — HIRES 60 fim. image (linear stretch from - 185 MJy/sr) of KR 140 with peaks corre- 
sponding to IRAS point sources indicated by the white crosses. Contours show 1420 MHz emission 
(contours at 5 - 10 K at 1 K intervals). The IRAS point source in the lower right corner of the 
image was not in the field of view of our SCUBA observations. 

Fig. 2. — Jiggle maps of the NW part of the partial dust shell surrounding KR 140 at 850 /im 
(top; linear stretch from —0.06 to 0.1 Jy/beam) and 450 /um (bottom; linear stretch from —0.2 to 
0.05 Jy/beam; smoothed to 850 um beam size). Contours of 12 //m emission (20 - 32 MJy/sr at 
2 MJy/sr intervals) show that the 850 fxm emission is coming from the same region as the bright 
infrared emission. The field of view is approximately 2' in each case centered on 0:2000 = 2'^19°^44®, 
<^2000 = 61°10'45". The hexagonal shape is due to the arrangement of the bolometers in the detector 
array. The location of these maps is indicated on the larger scan maps shown in Figure 3. 

Fig. 3. — Raw scan maps of KR 140 at 850 /xm (top; linear stretch from —0.083 to 0.151 Jy/beam) 
and 450 /Ltm (bottom; linear stretch from —0.62 to 0.86 Jy/beam). The locations of the jiggle maps 
shown in Figure 2 are indicated by the hexagons. 

Fig. 4. — Scan maps of KR 140, Weiner filtered and with long wavelength modes suppressed. Top 
shows 850 um emission (linear stretch from —0.045 to 0.11 Jy/beam), and bottom shows 450 /xm 
emission (linear stretch from —0.30 to 0.38 Jy/beam). Numbers and white crosses correspond to 
the source numbers and positions listed in Table 2. 

Fig. 5. — Scan map of KR 140 at 850 /xm (Figure 4) compared with 1420 MHz emission (top left, 
note central depression), HIRES 60 fim emission (top right), and CO (J=l-0) emission (bottom). 
1420 MHz contours are from 5 - 10 K at 1 K intervals and 60 ^m contours are from 100 - 250 MJy/sr 
at 25 MJy/sr intervals. The bottom left figure shows integrated CO emission from —45.5 km 
to —47.2 km s~^ (Vlsr); contours axe from 2.5 - 20 K at 2.5 K intervals. The bottom right figure 
shows integrated CO emission from —48.0 km to —53.7 km s~^ (Vlsr); note the doubling of 
scale with contours 5 - 40 K at 5 K intervals. 

Fig. 6. — Effect of changing temperature on the mass derived using equation 1. Deviations in 
the derived mass from starting temperatures of 10 K (solid line), 15 K (dashed line), and 25 K 
(dash-dot line) are shown. Horizontal solid lines indicate changes by a factor of two. At lower 
temperatures the derived mass is much more sensitive to temperature. 

Fig. 7. — SED for the sources KMJB 15, 3, and 1, constructed using cospatial HIRES and SCUBA 
data. The various symbols indicate measured fliix densities and the lines indicate the best-fitting 
bluebody curves (fit to submillimeter far-infrared points only). 

Fig. 8. — Closeup of the region surrounding IRAS 02168-1-6052. White crosses indicate KMJB 6 to 
8 (top to bottom). The star indicates the position of BIRS 128. The greyscale image is the 850 /xm 
scan map (—0.04 to 0.08 Jy/beam; linear stretch). Solid contours are HIRES 60 /um emission (130 
- 220 MJy/sr at 30 MJy/sr intervals). Dashed contours show integrated CO (J=l-0) emission 
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(contours at 7, 9 and 11 K; integrated from —45.5 km s ^ to —47.2 km s ^ (Vlsr))- 

Fig. 9. — Closeup of the region surrounding IRAS 02157+6053. White crosses indicate KMJB 17 
to 19 (1 - r). The greyscale image is the 850 /xm scan map hnearly stretched from —0.028 to 0.092 
Jy/beam. Contours of HIRES 60 /xm emission are overlayed (soUd contours, 130-180 MJy/sr at 10 

MJy/sr intervals). Integrated CO emission 

(Vlsr = -48.46 to -50.93 km s"^) is indicated by the 
dashed contours (6 - 7.5 K at 0.5 K intervals). 

Fig. 10. — Mass function for the submillimeter sources in KR 140. The cumulative clump number is 
plotted against mass (jagged solid line). Overlaid are power-law distributions with different indices, 
0.49 (best-fit; lower solid line), 1.0 (dash-dot line), and 1.5 (steeper like typical stellar IMF; dash 
line). The 1.0 and 1.5 curves have been fixed to match the best-fit curve at the third most massive 
object. 

Fig. 11. — Like Figure 10 but for the large cores in p Oph (Motte et al. 1998). Distributions with 
power-law indices 0.60 (best-fit; solid lower line), 1.0 (dash-dot) and 1.5 (dash) are shown. 

Fig. 12. — Mass determined from the submillimeter flux density versus the Bonnor-Ebert mass 
(Mbe) calculated for a sphere of the same radius and temperature. The solid diagonal line (line 
of instability) represents the condition M = Mbe- Crosses indicate sources with M > Mbe and 
diamonds indicate those with M < Mbe- The asterisks represent the two sources known to contain 
embedded stars. Numbers correspond to the sources shown in Figure 4. 

Fig. 13.— Column density M/D^^ versus linear size. The equivalent value of total visual extinction 
is shown on the right y-axis. Note the bimodal distribution of sources. S3mibols are the same as in 
Figure 12. 
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Table 1. IRAS Point Sources Near KR 140 
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Table 2. Submillimeter Sources 
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'^Acronym for the new submillimeter sources; from Kerton, Martin, Johnstone, &: Ballantyne 

from submm-IR SED fit 
'^T assumed 

from 450/850 /xm color 
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Table 3. Dust Opacity at 850 fxm 



C (g cm"2)* 


Source 


Recommended Use 


286 


Draine (1990) & Draine & Lee (1984) 


diffuse ISM 


116 


Hildebrand (1983) 


diffuse ISM, canonical value 


85.5 


Preibisch et al. (1993) 


dense clumps and cores 


53.4 


Ossenkopf &l Henning (1994) 


regions with n ~ 10^ cm^ 


38.9 


Ossenkopf k, Henning (1994) 


regions with n > 10^ cm^ 


21.4 


Kriigel & Siebenmorgen (1994) 


dense, cold regions 



combines the dust opacity {kv) and the dust to gas mass ratio (C = Mg/^M^Ki,)) 
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